Stabilization of a (3+l)D soliton in a Kerr medium by a rapidly oscillating dispersion 
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Using the numerical solution of the nonlinear Schrodinger equation and a variational method 
it is shown that (3+l)-dimensional spatiotemporal optical solitons can be stabilized by a rapidly 
oscillating dispersion coefficient in a Kerr medium with cubic nonlinearity. This has immediate 
consequence in generating dispersion-managed robust optical soliton in communication as well as 
possible stabilized Bose-Einstein condensates in periodic optical-lattice potential via an effective- 
mass formulation. We also critically compare the present stabilization with that obtained by a rapid 
sinusoidal oscillation of the Kerr nonlinearity parameter. 
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I. INTRODUCTION 

After the prediction of self-trapping pj of an optical 
beam in a nonlinear medium resulting in an optical soli- 
ton 0, 0] , there have been many theoretical and exper- 
imental studies to stabilize such a soliton under differ- 
ent conditions of nonlinearity and dispersion. A bright 
soliton in (1+1) dimension (D) in Kerr medium (cubic 
nonlinearity) is unconditionally stable for positive or self- 
focusing (SF) nonlinearity in the nonlinear Schrodinger 
equation (NLS) 0. However, in (2+l)D and (3+l)D, 
in homogeneous bulk Kerr medium one cannot have an 
unconditionally stable soliton-like beam 0] . (If the non- 
linearity is negative or self-defocusing (SDF), any ini- 
tially created soliton always spreads out Q.) If the non- 
linearity is positive or of SF type, any initially created 
highly localized soliton is unstable in both (2+l)D and 
(3+l)D 3]. Such a confined wave packet in (3+l)D is 
often called a light bullet and represents the extension 
of a self-trapped optical beam into the temporal domain 
• The generation of a light bullet is of vital importance 
in soliton-based communication systems. On the experi- 
mental front strong stabilization of (2+l)D discrete vor- 
tex solitons in a SF nonlinear medium have been observed 
in optically induced photonic lattices Q|. Such solitons, 
however, can only be modeled by a more complicated 
nonlinearity in the NLS in contrast to the simple cubic 
Kerr nonlinearity considered in this paper. 

Recently, through a numerical simulation as well as 
a variational calculation based on the NLS it has been 
shown that a soliton in (2+l)D 0,0] can be stabilized in 
a layered medium if a variable Kerr nonlinearity coeffi- 
cient is used in different layers 0, @ ■ A weak modulation 
of the nonlinearity coefficient along the propagation di- 
rection leads to a reasonable stabilization in (2+l)D 5]. 
A much better stabilization results if the Kerr coefficient 
is a layered medium is allowed to vary rapidly between 
successive SDF and SF type nonlinearities, i.e., between 
positive and negative values 0. 

However, in recent years optical soliton formed by a 
controlled variation of dispersion coefficient (dispersion- 
managed optical soliton) in a Kerr medium has been 
of general interest in communication 7J. A dispersion- 
managed soliton allows robust propagation of pulses and 



is favored over normal solitons. The recent stabilized op- 
tical beam in a nonlinear waveguide array in (1+1)D [||, 
called diffraction managed soliton, is quite similar to the 
present dispersion-managed optical soliton. That prob- 
lem was described |8( by the one dimensional NLS with a 
periodic dispersion coefficient. Such a periodic variation 
of the dispersion coefficient leads to a greater stability 
of the soliton during propagation in (1+1 )D Q. More 
recently, Abdullaev et al Q have shown by a variational 
and numerical solution of the NLS that by employing a 
rapidly varying dispersion coefficient it is possible to sta- 
bilize an optical soliton in (2+l)D over large propagation 
distances. 

The possibility of the stabilization of an optical soliton 
in (2+l)D by dispersion ma nag ement [t| and Kerr sin- 
gularity management [lfllTlllial as well as in (3+l)D by 
Kerr singularity management |13l Il4| in SF medium has 
led us to investigate the possibility of the stabilization of 
an optical soliton in (3+l)D by dispersion management. 
However, some of these studies were concerned with the 
stabilization of Bose-Einstein condensates us- 
ing the nonlinear Gross-Pitaevskii equation |15| which is 
very similar to the NLS structurally. The extension of 
the problem of stabilization of optical soliton from lower 
to higher dimensions, e.g. from (1+1)D to (3+l)D, is of 
great interest in the actual three-dimensional world. It 
is also more challenging due to the possibility of violent 
collapse in SF medium in higher dimensions 

We begin the present study with a time-dependent 
variational solution of the NLS with a rapidly oscillating 
dispersion coefficient in (3+l)D. The variational method 
leads to a set of coupled differential equations for the soli- 
ton width and chirp parameters, which is solved numeri- 
cally by the fourth-order Adams-Bashforth method [16j. 
The variational solution illustrates that a (3+l)D light 
bullet can be stabilized over a large propagation distance 
in a Kerr SF medium with a rapidly oscillating dispersion 
coefficient for beam power above a critical value. Next we 
turn to a complete numerical solution of the partial dif- 
ferential NLS by the Crank-Nicholson method [13 ■ The 
complete numerical solution demonstrates the stabiliza- 
tion of a light bullet over large propagation distance by 
dispersion management in (3+l)D SF Kerr medium. 

We also compare the present stabilization of a light 
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bullet by dispersion management with the stabilization 
obtained by applying an oscillating Kerr nonlinearity as 
suggested in Ref. Q. We find that both schemes leads 
to comparable stabilization of a light bullet over a large 
propagation distance of few thousand units. 

Although, the present work is of primary interest in 
the generation of robust optical solitons, it is also of 
interest from a mathematical point of view in nonlin- 
ear physics in the stabilization of a soliton in three di- 
mensions. Moreover, this investigation has interesting 
implication in the study of Bose-Einstein condensates. 
The quantum nonlinear equation governing the evolu- 
tion of the condensate, known as the Gross-Pitaevskii 
equation [l^ , is identical with the (classical) NLS Q , for 
the evolution of optical soliton, although the interpreta- 
tion of the different variables of these two equations is 
different. In recent years there has been experimental 
[l8l | and theoretical pij studies of Bose-Einstein conden- 
sates in periodic optical-lattice potentials generated by a 
standing-wave laser beam. By an effective-mass descrip- 
tion, the Gross-Pitaevskii equation for the condensate in 
a periodic optical-lattice potential can be reduced to a 
dispersion-managed NLSjWhere the effective mass could 
be positive or negative 20]. There has already been ex- 
perimental application of such dispersion management to 
Bose-Einstein condensates |2^. The effective mass can 
be varied by changing the parameters of the periodic po- 
tential 9]. Once a controlled variation of the effective 
mass would be possible, it might be possible to stabilize 
a Bose-Einstein condensate in this fashion. 

In Sec. II we present a variational study of the prob- 
lem and demonstrate the possible stabilization of a light 
bullet in (3+l)D. In Sec. Ill we present a complete nu- 
merical study based on the NLS to study in detail the 
stabilization of a dispersion-managed light bullet. We 
also compare the present stabilization with that by non- 
linearity management as in Ref. [14(. Finally, in Sec. IV 
we give the concluding remarks. 
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II. VARIATIONAL CALCULATION 



For anomalous dispersion, the NLS can be written as 



;9_ 7(2), 

'dz + 2 



\u(r,z)\< 



u(r,z) = 0, 



(1) 



where in (3+l)D the three dimensional vector r = 
(x,y,t) has space components x and y and time com- 
ponent t, and z is the direction of propagation. The 
Laplacian operator "S7 2 acts on the variables x, y, and 
t. The dispersion coefficient 7(2) in a Kerr medium is 
taken to be rapidly oscillating and can have successive 
positive and negative values. The normalization condi- 
tion is J dr\u(r, z)\ 2 = P, where P is the power of the 
optical beam 0, |(| . 

For a spherically symmetric soliton in (3+l)D, 
u(r,z) = U(r,z). Then the radial part of the NLS 



becomes [J 
d 7 (» d 2 



dz 



2 dr 2 



700 d 
r dr 



\U{r,z)\' 



U(r,z) = Q. (2) 



For stabilization of the soliton, we shall employ the fol- 
lowing oscillating dispersion coefficient: 7(2) = [1 + 
gi sin(wz)], where frequency u> is large. In this paper we 
consider variational and numerical solutions of Eq. © 
to illustrate the stabilization of a (3+l)D light bullet. 

Recently, it was suggested jlj| that a rapid variation of 
the Kerr nonlinearity parameter also stabilizes a soliton. 
In this study we critically compare these two ways of 
stabilization. For that purpose we consider 
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U(r,z) = 0, (3) 



where the nonlinearity parameter will be taken to be 
rapidly varying for stabilizing the soliton. In Ref. [l4j 
n(z) was taken to be rapidly jumping between positive 
and negative values and a small propagation distance z 
was considered for stabilization. In the present study 
a sinusoidal variation of n(z) is considered over a much 
larger propagation distance. For stabilization we consider 
n(z) — [1 — <72 sin(wz)] in this paper, where the second 
term on the right is rapidly oscillating for a large u> and 
where gi is a constant. 

First we consider the variational approach with the 
followingtrial Gaussian wave function for the solution of 
Eq. © BII3IH0 



U(r, z) = N(z) exp 



2R 2 {z) 



-b(z)r 2 



- ia(z) 



, (4) 



where N(z), R(z), b(z), and a(z) are the soliton's am- 
plitude, width, chirp, and phase, respectively. In Eq. Q 
N(z) = P^ 2 /[tt 3 / 4 R 3 / 2 (z)}. The trial function © satis- 
fies (a) the normalization condition 0, 0] as well as the 
boundary conditions (b) U(r, z) — > constant as r — > 
and (c) \U(r, z)\ decays exponentially as r ooJQ. 
The Lagrangian density for generating Eq. @ is 0, 

M 



C(U) = 



dU^ 

d; 



-U* 



9U1 U )_1 

dz 2 



dU 



\u\ 4 



(5) 



where the functional dependence of the quantities on r 
and z is suppressed. The trial function Q is substi- 
tuted in the Lagrangian density and the effective La- 
grangian is calculated by integrating the Lagrangian den- 
sity: L eS = j C(U)dr. With U given by Eq. Q the 
effective Lagrangian is 

v(f = \^ /2 N 2 {z)R\z) 
3 j(z) 



L, 



2 

N 2 (z) 

2^2 ~ 2^ 



- lkz)R 2 (z) 



2a(z) 



-R 2 (z)b 2 (zh(z) 



(6) 



where the overhead dot denotes derivative with respect 
to z. The standard Euler-Lagrange equations for N(z), 
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R(z), b(z), and a(z) are HIH0 given by 

d d£ e ff = dL eS 

dz dp 9/3 ' 



(7) 



where (3 stands for N(z), R(z), b(z), or a(z). The Euler- 
Lagrange equation for a(z) gives the normalization: 



n 3/2 N 2 R 3 = P. 



(8) 



The Euler-Lagrange equations for R(z) and N(z) are 
given by 



ri 4 " 7 ,,2 ^ 2 



(9) 
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where the 2 dependence of the variables is suppressed. 
The Euler-Lagrange equation for b{z) is 



R(z) = R(z)b(z)~f{z). 



(11) 



Equations © and 11U|) leads to the following rate of vari- 
ation of b{z): 



Kz) = ^ ) -b 2 (z) 1 {z) 



(12) 



Now defining the new variable k(z) = b(z)R(z), Eqs. 
(|ll|l and l|12[l can be rewritten as 



R(z) 
k(z) 



j(z)k(z) 
l(z) P 

R 3 (z) 2VYi 



■R\zY 



(13) 
(14) 



Equations (|13fl and (|14|l are the variational equations of 
motion for R(z) and k(z). Here we consider the stabil- 
ity condition of optical solitons described variationally 
by these equations for 7(2) = [1 + gi sm(ujz)], where the 
second term is rapidly varying (large lo) with zero mean 
value. However, this set of coupled differential equa- 
tions cannot be solved analytically and we resort to their 
numerical solution using the four-step Adams-Bashforth 
method ^(|. For a large u> and P, and g\ > 2 stabiliza- 
tion of the soliton is always possible provided that the 
initial values of R(z) and k(z) are appropriately chosen. 
The soliton can also be stabilized by employing a peri- 
odically oscillating Kerr nonlinearity. However, in that 
case the coupled set of differential equations for width 
R{z) and chirp b(z) was simpler and could be studied 
analytically to establish stabilization of the system [T^j . 

In Fig. 1 we exhibit the numerical results for the width 
R(z) obtained from the solution of the variational equa- 
tions l|13f) and (|14fl . In the present calculation we take 
to = 15, gi = 4, P = 500, and R(z = 0) = 6 If the ini- 
tial value of b(z) is chosen appropriately, stabilization of 
the system can be obtained. For these set of parameters, 
stabilization was obtained for b(z = 0) = 0.14, and we 



FIG. 1: The width of the soliton R(z)/R(0) vs. propagation 
distance z from the numerical solution of the coupled set of 
equations 11311 and 1141 . with 7(2) = [1 + gi sin^z)] for P — 
500, gi =4,cj= 15, R(0) = 6 and k(0) = 0.14. 



show the result in Fig. 1, where we plot R(z)/R(0) vs. z. 
The width R{z) is found to to remain stable for a range 
of values of z. Small oscillations in Fig. 1 sustain the 
soliton for a large propagation distance z. The numeri- 
cal values for g\ and lo are taken as examples, otherwise 
they do not have great consequence on the result so long 
as to is large corresponding to rapid oscillation. There is 
no upper limit for P and stabilization seems possible for 
an arbitrarily large P. For larger P the stabilization is 
more sustained and consequently, it is easier to stabilize 
a soliton. However, SF nonlinearity is essential for the 
stabilization and stabilization is only possible for power 
P above a critical value. 

The sinusoidal variation of the dispersion coefficient as 
considered above in the variational study only simplifies 
the algebra and is by no means necessary for stabiliza- 
tion of solitons. Any rapid periodic oscillation of the 
dispersion coefficient was also found to stabilize the soli- 
ton. The same was found to be true in the stabilization 
of a soliton by Kerr-singularity management in (3+l)D 
00. 



III. NUMERICAL RESULT 



Motivated by the approximate variational study above, 
now we solve Eqs. @ and © numerically using a split- 
step iteration method employing the Crank-Nicholson 
discretization scheme |0, |23| • The ^-iteration is started 
with the known solution of some auxiliary equation with 
zero nonlinearity. The auxiliary equation with known 
Gaussian solution is obtained by adding a harmonic os- 
cillator potential r 2 to Eqs. © and J3J and setting the 
nonlinear term to zero, e.g., 
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d(z)r 



U(r,z)=0. (15) 
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FIG. 2: (Color online) Variation of power P(z) and strength 
d(z) of radial trap in the initial stage of stabilization until the 
radial trap removed and the final power P — 868.4 attained 
at 2 = 400. 



In Eq. Ijl 5(1 we have introduced a strength parameter 
d(z) with the radial trap. Normally, d(z) = 1; when 
the radial trap is switched off d(z) will be reduced to 
zero. Asymptotically decaying real boundary condition 
for U(r, z) is used throughout this calculation. 

When the radial trap is reduced to zero the wave func- 
tion extends to a large value of r. Hence to obtain a con- 
verged result a large value of r is to be considered in the 
Crank-Nicholson discretization scheme. In the present 
calculation we employ an r-step of 0.1 with 40 000 grid 
points spanning an r-value up to 4000 and a z-step of 
0.001. In the course of z-iteration a positive SF Kerr 
nonlinearity with appropriate power P is switched on 
slowly and the harmonic trap is also switched off slowly. 
If the nonlinearity is increased rapidly the system col- 
lapses. The tendency to collapse or expand to infinity 
must be avoided to obtain a stabilized soliton. Although, 
for the sake of convenience we applied a harmonic trap in 
the beginning of our simulation, which is removed later 
with the increase of nonlinearity, this restriction is by no 
means necessary for stabilizing a soliton. 

First we consider the numerical solution of Eq. (J2J) to 
create the desired soliton. After switching off the har- 
monic trap and introducing the final power P, the os- 
cillating dispersion coefficient j(z) = [1 + g\ sin(cjz)] is 
introduced. A stabilization of the final solution could 
be obtained for a large u> and P. If the SF power af- 
ter switching off the harmonic trap is large compared to 
the spatiotemporal size of the beam the system becomes 
highly attractive in the final stage and it eventually col- 
lapses. If the final power after switching off the harmonic 
trap is small for its size the system becomes weakly at- 
tractive in the final stage and it expands to infinity. The 
final power has to have an appropriate value consistent 
with the size, for final stabilization. 

Starting from Eq. ((15(1 the soliton can be created with 
a wide range of variation of power P(z) and strength d{z). 
In Fig. 2 we show the actual variation of P(z) and d(z) 




Dispersion Management 




Nonlinearity Management 




FIG. 3: (Color online) (a) The solitonic wave function 
\U(r, z)\ vs. r and z of a freely propagating light bullet of 
power P = 868.4, the arrow at z = 300 shows the begin- 
ning of collapse with \U(r — 0,z)\ starting to diverge; (b) 
The solitonic wave function \U(r, z)\ of Eq. (|5J vs. r and 
z of the light bullet of (a) stabilized by dispersion manage- 
ment 7(2;) — [1 + 4sin(57rz)] applied for all z > 0; (c) The 
solitonic wave function \U(r, z)\ of Eq. © vs. r and z of the 
light bullet of (a) stabilized by Kerr nonlinearity management 
n(z) = [1 — 4sin(57rz)] applied for all z > 0. 
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of the radial trap in Eq. (|15fl . A stabilized soliton can be 
obtained for different variations of P(z) and d(z) and the 
variation in Fig. 2 are by no means unique and are to be 
considered as a possible variation. By switching off the 
harmonic trap linearly for 13 < z < 15 and increasing 
the nonlinearity for < z < 400, an almost stationary 
light bullet is prepared. One interesting feature of Fig. 2 
is that we could not find a linear variation of P(z) over 
the whole range to obtain the final soliton. A linear vari- 
ation usually led to collapse or expansion to infinity of 
the final soliton. A fine-tuning of the power was needed 
for obtaining the quasi-stationary soliton with the kinetic 
pressure almost balancing the nonlinear attraction. This 
fine-tuning led to the final fractional power P = 868.4. 
A quasi-stationary soliton, rather than a rapidly expand- 
ing or collapsing soliton, was ideal for stabilization. Al- 
though the final power P = 868.4 is used in this study, 
stabilization can be obtained for any P greater than a 
critical value. 

In the preparation of the soliton in Fig. 2 no oscillat- 
ing dispersion coefficient has been used. The oscillating 
dispersion coefficient has been applied later. In (2+l)D 
Saito and Ueda first applied a weak oscillating non- 
linearity and then increased its strength in a linear fash- 
ion with z iteration. The nonlinearity manipulation of 
Saito and Ueda ^2] is distinct from ours. However, the 
final result should be independent of how the oscillating 
dispersion or nonlinearity coefficient is introduced. 

The propagation of the free light bullet is studied by 
a direct numerical solution of the NLS J2J- Once it is 
allowed to propagate in the z direction the light bullet 
starts to collapse and is destroyed after propagating a 
distance of about 300 units of z as shown in Fig. 3 (a) 
where we plot the profile of the freely propagating light 
bullet |J7(r, z)\ vs. r and z. This means that the non- 
linear attraction supersedes the kinetic pressure in this 
soliton and it collapses eventually. 

In the present scheme of stabilization the dispersion co- 
efficient "f(z) = [l+^i sin(wz)] oscillates rapidly between 
positive and negative values. A negative dispersion coef- 
ficient 7(2) corresponds to collapse and a positive disper- 
sion coefficient corresponds to expansion. In one half of 
the oscillating cycle the soliton tends to collapse and in 
the other half it tries to expand. If, in the first semi-cycle 
of the application of the oscillating dispersion coefficient, 
7(2:) is positive, the slowly collapsing soliton of Fig. 3 
(a) will tend to expand in this semi-cycle. In the next 
semi-cycle, 7(2) will be negative and the soliton will tend 
to collapse. This happens for a positive gi > 1. If the 
expansion in the first interval is compensated for by the 
collapse in the next interval, a stabilization of the system 
is obtained. However, if gi is negative, in the first semi- 
cycle of oscillation the slowly collapsing soliton of Fig. 3 
(a) will collapse further and a satisfactory stabilization 
cannot be obtained. For a large u the system remains vir- 
tually static and the very small oscillations arising from 
collapse and expansion remain unperceptible. 

Next using Eq. © we illustrate the possible stabiliza- 
tion of the (3+l)D soliton of power P = 868.4 of Fig. 3 
(a) by the application of the oscillating dispersion coeffi- 
cient 7(2) = [1 + .gi sin(wz)] with g\ — 4 and uj — 57T for 
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FIG. 4: (Color online) The wave function \U{r,z = 1000) | of 
the stabilized light bullet of Figs. 3 (b) and (c) vs. r. For 
comparison the Gaussian (~ exp(— r 2 /100)) is also shown. 



z > 0. The profile of the dispersion-managed light bullet 
\U(r, z)\ vs. r and z is plotted in Fig. 3 (b). It is realized 
that after the application of dispersion management the 
soliton remains stable for a large propagation distance z. 
This shows clearly the effect of the oscillating dispersion 
coefficient on stabilization. The dispersion management 
of the type considered in Fig. 3 (b) can prolong the life 
of the soliton significantly. In Fig. 3 (b) the profile of 
the solitonic wave function over the large interval of z 
clearly shows the quality of stabilization. The stabiliza- 
tion seems to be good and can be continued for longer 
intervals of z by considering a soliton of larger power. 

Finally, using Eq. 10 we consider the stabilization of 
the above soliton of power P = 868.4 by the application 
of an oscillating Kerr nonlinearity n{z) = [1 — 32 sin(u;z)] 
with gi = 4 and lj = 5n for z > 0. In this case for 
gi > 1, the Kerr nonlinearity coefficient n(z) is negative 
in the first semi-cycle of the oscillating Kerr coefficient. 
This will correspond to an expansion of the soliton, so 
that the collapsing tendency of the soliton of Fig. 3 (a) 
could be stopped and a stabilization of the soliton could 
be obtained. The profile of the stabilized light bullet by 
nonlinearity management is plotted in Fig. 3 (c). It is 
found that the oscillating Kerr nonlinearity can stabilize 
the light bullet for a large propagation distance and en- 
hance the life of the soliton significantly. The quality of 
stabilization is comparable to that obtained by disper- 
sion management. In both schemes the light bullet can 
remain stable for few thousands units of z. In Figs. 3 we 
illustrate this stabilization for z up to 2000. The stabi- 
lization obtained i n (2+l) D by similar dispersion and 
nonlinearity [61 flfl fill Il2| managements was over a much 
smaller interval of z (< 100). 

The final solitonic wave function \U(r, z)\ extends to 
very large values of r. It has a sharp peak for small r, 
reminiscent of strong binding, and a large extention in 
r, typical of weak binding. We illustrate these features 
in Fig. 4 where we plot \U(r, z = 1000) | of a light bul- 
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let vs. r for dispersion or nonlinearity managements of 
Figs. 3 (b) and (c). The curves for dispersion and non- 
linearity managements are identical with each other and 
a single curve is shown. The long tail of the stabilized 
light bullet is shown in inset. For comparison a Gaussian 
(~ exp(— r 2 /100)) is also shown. 

For both dispersion and nonlinearity managements we 
found that the stabilization can only be obtained for 
beams with power larger than a critical value. However, 
we could not obtain a precise value of this critical power 
numerically. In Ref. [14| we obtained a variational esti- 
mate of this critical power (~ 40) for nonlinearity man- 
agement. We could not obtain such an estimate for dis- 
persion management. Numerically, we found that it was 
easier to obtain stabilization of beams with power much 
larger than the critical value. 



IV. CONCLUSION 

In conclusion, after a variational and numerical study 
of the NLS we find that it is possible to stabilize a spa- 



tiotemporal light bullet by employing a rapidly oscillat- 
ing dispersion coefficient. We find that the nature of 
this stabilization is similar to that obtained by a rapid 
variation of the Kerr-nonlinearity parameter. In both 
cases stabilization is possible in a SF Kerr medium for 
optical beams with power P larger than a critical value. 
However, there is no upper limit of P for stabilization. 
A larger value of P is preferred for stabilization. The 
stabilization of the light bullet by dispersion and Kerr 
nonsingularity management ,14| should find experimen- 
tal application in optics and Bose-Einstcin condensation. 
Such a dispersion- and nonsingularity-managed optical 
soliton can propagate large distances with a minimum of 
distortion and is to be preferred over normal solitons in 
optical communication. 
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